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Abstract 

In the context of the sensational results concerning superluminal velocities, 
announced recently by the OPERA Collaboration, we have proposed a clas- 
sical model yielding a statistically calculated measured velocity of a beam, 
higher than the velocity of the particles constituting the beam. The two key 
elements of our model, necessary and sufficient to obtain this curious result, 
are a time-dependent "transmission" function and statistical method of the 
maximum-likelihood estimation. 
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Inspired by the amazing results of the OPERA Collaboration claiming 
that the velocity of light c has been beaten by a beam of neutrinos, we pro- 
pose a simple classical model, which yields the curious effect of a seeming 
increase of "effective" velocity (see for an earlier independent approach, 
and also [3], for a wave version). There are dozens of papers, which have 
recently appeared, adopting various attitudes towards the results presented 
by the OPERA Collaboration, for example [H], [g], 0, [sl (see also an example 
of an older work on superluminal velocities j9|). The attitude of our paper is 
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definitely skeptical. In view of the fact that we have found a ("mathematical") 
model providing an "artificial" increase of real velocity, we are forced to put 
the results announced by the OPERA Collaboration in doubt. Moreover, we 
have also proposed a simple working example, such that appropriately fit- 
ting its parameters one can simulate the controversial results of the OPERA 
Collaboration. Strictly speaking, our paper does not indisputably invalidates 
the conclusions drawn by the OPERA Collaboration, but it seriously weakens 
their argumentation, indicating a logical gap in their reasoning. 

Our model is purely classical and dynamics free. No new physics, nor quan- 
tum mechanics, nor even (classical) wave mechanics is involved. Only stan- 
dard classical kinematics notions as well as the statistical method of the 
maximum-likelihood estimation (MLE) are used in our approach. 

The assumptions of our model are the following. A spatially homogeneous, 
lasting the period T, beam of classical particles ("extraction") moving with 
a constant speed -u (^ c) travels from a source to a detector. The distance 
between the source and the detector is d. The probability density function 
(PDF) of the time of emission of the particles within the duration T of pro- 
duction of the beam is given by the function w{t). In an ideal situation (none 
of the particles is lost) we would obtain, according to classical kinematics, the 
measured data waveform y{t) = w{t — to), where to = d/u. Now, let as sup- 
pose that the fraction of the particles emitted, measured by the detector, due 
to some physical mechanism, is given by the (non-negative) "transmission" 
function 

where (t) is the number of the particles emitted at the time instant t, and 
Nd (t) is the number of the particles detected, which have been sent at the 
same time instant t. In other words, in general, not all particles emitted are 
detected (fit) < 1) — obvious, and moreover f(t) 7^ const — conceivable. 
For simplicity, we will assume that the numbers of the particles, N(i{t) and 
iVe(t), are so large that we are allowed to use a continuous approximation. 
Then, the transmission function f{t) is a (continuous) function satisfying the 
condition ^ /(t) ^ 1. 

To be able to draw conclusions from experimental data, we should implement 
some statistical methodology. To be so precise as possible, in our approach, 
we adopt the method of the MLE, as the OPERA Collaboration has done 
hj. In the framework of the MLE, we introduce the likelihood function L. 
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The logarithm / of L is given by the formula 



m ^ log [im ^ Yl [^(^i + '^^)] ' (2) 

where tj are the time instants corresponding to the measurement events at 
the detector, and the time deviation 6t we are interested in (see [1], for de- 
tails) provides the maximum of / (and also of L). As the numbers of the 
measured events tj are large, in our continuous approach, instead of summa- 
tion in ([2]), we should use integration with an appropriate integration measure 
z{t)dt, where z{t) represents the time distribution of the experimental events 
detected by the OPERA. In fact, z{t) is determined by product of two fac- 
tors. The first factor, y, is proportional to the number of the particles sent, 
i.e. y{t + to) = w{t), and the second one is proportional to the transmission 
function f{t). Then, 

zit) = f{t)y{t + to) = f{t)w{t). (3) 

Finally, 

l(6t) = J log [w{t + 6t)] f{t)w{t)dt. (4) 

To demonstrate that our idea actually works, we propose a specific example. 
The parameters of the example are so fitted that it yields the time deviation 
5t ^ +75.5 ns (for comparison, the OPERA result is 60.7ns). For calcula- 
tional simplicity, we assume the Gaussian form of the PDF (see the solid line 
in FigiH), 

w{t) = exp ' (5) 

as well as the Gaussian form of of the transmission function, 

/W = l-^exp[-(t-l)^], (6) 

where the one tenth in front of the exponent is reminiscent of "10% variation" 
in The time distribution of the "detected experimental data" z{t), corre- 
sponding to w{t) of the form ([5]) and to the transmission function f{t) of the 
form ([6]) is presented in Fig{T]by the dashed curve. In this (doubly) Gaussian 
case it is even possible to solve the problem analytically (see Appendix), but 
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Figure 1: 

The presented specific example is intended to "mathematically" qualitatively 
simulate the OPERA experiment, yielding the time deviation 6t ~ +75.5 ns. The 
solid curve represents the probability density function (PDF) w{t), whereas the 
dashed one corresponds to the time distribution of the "detected experimental 

data" z{t). 

for our purposes a numerical value will do. It appears, that the maximum of 
(jl]) is attained for the time deviation 

5t ^ +0.0288 s.d., (7) 

where s.d. means the standard deviation. One can easily translate the 
dimensionless ([7]) into a dimensionfull entity. In the OPERA experiment 
T ?a 10.5 /is, whereas in our example we can reasonably assume, for definite- 
ness, that T equals (in dimensionless units, or in the units of the standard 
deviation) twice the two standard deviations, i.e. T ^ 2 ■ 2 = 4 (see two verti- 
cal intervals in Fig{T]). Then, the dimensionfull time deviation corresponding 
to dZD is 

0X)288^KL5^s^ 75.5 ns ^ 60.7ns (OPERA). (8) 

In particular, our analysis confirms the findings of Q that the time devia- 
tion 6t is independent of the distance d, but depends on the shape of the 
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beam. Obviously, the numerical coincidence (jSj) has only an illustrative pur- 
pose. The non-zero time deviation 5t directly implies a modified velocity v 
according to the elementary formula 



In our model the sign of 5t (positive in our example) depends on details of 
the form of the transmission function f{t). E.g., reversing the sign in front 
of the exponent in ([6]) reverses the sign of 5t. 

We would like to stress that our specific example is not supposed to mimic a 
real situation in the OPERA experiment, but only to demonstrate that it is, 
in principle, possible to easily come to its conclusions without any reference to 
superluminal particles and/or some exotic phenomena. Any considerations 
concerning a possible physical mechanism governing the time-dependence of 
the "transmission" function f{t) are outside the scope of our paper. We can 
only speculate that f{t) ^ const could be a property of the source, or of the 
detector or it could follow from interactions in the Earth's crust. Moreover, 
measurement errors, always encountered in real experiments, do not enter 
our considerations as they have nothing to do with the discussed effect. 

It is also possible to intuitively explain the non-zero value of 5t. Namely, de- 
forming the shape of the PDF w{t) with an appropriate (time asymmetric) 
transmission function f{t) shifts a portion (its "upper" part) of the "wave- 
form" w{t) forward or backward, yielding positive or negative 5t, respectively, 
which is next erroneously interpreted by the method of the MLE as a "real" 
time deviation. In a sense, our approach is a "corpuscular" analog of the 
well-known situation of "superluminal" velocity of light in material media 
flo| . where to some extent, the role of f{t) plays dispersion. 

In conclusion, we would like to emphasize that there are two elements of 
our model, responsible for the curious effect of the superluminal velocity: 
the time-dependent transmission function /(t), which should "favor" leading 
edges and "disfavor" trailing edges of the beam, and statistical methodology, 
erroneously assuming the method of the MLE. This is why it is conceivable 
that the paradoxical results of the OPERA Collaboration could be possibly 
avoided upon another approach, e.g. an approach giving preference to the 
front velocity or to another statistical method. Thus, no superluminal par- 
ticles are necessary to explain the "superluminal" velocities derived by the 
OPERA Collaboration. Moreover, independently of the further fate of the 
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conclusions drawn by the OPERA Collaboration (acceptance or refutation), 
we have demonstrated that standard statistical methods, e.g. MLE, should 
be used with great they can provide erroneous output. Therefore, the 

results of our paper do not rely on the final solution of the OPERA paradox. 



Appendix A. Analytical form 

For the Reader's convenience, we present here analytical form of the expres- 
sions used in the main part of the paper. For w{t) given by ([5]) and f{t) 
given by (jS]), we obtain the integral (jl]) 
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exp [-j]dt 



_ y/^[V3{7 +12 6t + 9 6t^) - 270{/e{l + 6t' 
~ 270^2^ 

The maximum of l{6t) is determined by its derivative 

[VS (12 + 18 St) - bAO^St] 



I' {6t) 
and is attained for 



270^2^ 



(A.l) 



(A.2) 



6t 



3073^-3' 



(A.3) 
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